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1.1 Introduction 

One of the main goals of Physics is to identify and characterize the possi- 
ble phases of matter, as well as the transitions between these phases. In 
this chapter we are concerned with theoretical and computational aspects of 
quantum phases and quantum phase transitions, involving extended quantum 
many-body systems at zero temperature. 

Specifically, we consider the problem of constructing a real-space coarse- 
graining transformation for quantum lattice systems at zero temperature. One 
such transformation should dispose of the degrees of freedom corresponding 
to small distances while preserving the low energy properties of the system. 
Obtaining an effective description of a system in terms of less degrees of 
freedom is obviously very convenient for numerical calculations. Our priority 
here is, however, to explore the use of the coarse-graining transformation 
within the context of the renormalization group (RG) [1, 2, 3, 4]. Accordingly, 
we aim to define a RG flow in the space of local Hamiltonians, and to study 
quantum phases and quantum phase transitions by characterizing the fixed 
points of this flow. 

The coarse-graining scheme that we will describe is known as entanglement 
renormalization. Entanglement renormalization has been shown to be suitable 
to address the emergence of different types of order in systems of quantum 
spins in one and two spatial dimensions, including symmetry-breaking order 
and topological order, as well as to characterize quantum critical points. In 
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addition, the formalism can be generalized to study models where the basic 
degrees of freedom are fermionic or, more generally, anyonic. We will introduce 
the approach step by step throughout the chapter, and will then apply it to 
the characterization of quantum critical points. 

The content is organized in sections as follows. Sect. 1.2 discusses a coarse- 
graining transformation, based on isometries, that is capable of preserving 
ground state properties, but that fails to remove some of the short-distance 
degrees of freedom, which accumulate over successive iterations. This degrees 
of freedom are associated to short-range entanglement in the ground state 
of the system. Failure to properly get rid of short-range entanglement has 
important numerical implications, and precludes the use of RG ideas within 
this first coarse-graining scheme. 

Sect. 1.3 explains how to remove short-range entanglement from the lattice. 
This is achieved by introducing disentanglers that act across the boundary of 
blocks of sites before the coarse-graining step. In this way all short-distance 
degrees of freedom are removed from the system. Local operators are seen 
to be mapped into local operators by the ascending superoperator A, while 
its dual, the descending superoperator T>, allows us to recover a fine-grained 
reduced density matrix from a coarse-grained one. We also describe the class 
of states that can be exactly represented using the entanglement renormaliza- 
tion scheme, known as the multi-scale entanglement renormalization ansatz 
(MERA). 

Sect. 1.4 connects the present approach with the RG formalism. A RG 
map lZ(h) in the space of two-site interactions h is defined in terms of the 
ascending superoperator A. Each fixed point of the resulting (discrete) RG 
flow corresponds to a ground state that is invariant under the coarse-graining 
transformation and is described by a scale invariant MERA. The entangle- 
ment of scale invariant ground states allows to distinguish between different 
types of fixed points. At the fixed point, the ascending superoperator, now 
independent of the length scale, is known as the scaling superoperator S. 

Sect. 1.5 applies entanglement renormalization to the study of quantum 
phase transitions. Many universal properties of a quantum critical point, in- 
cluding critical exponents and the conformal data characterizing the pertinent 
conformal field theory, can be extracted from the scaling superoperator S. We 
also address boundary critical phenomena by adding a boundary to the scale 
invariant MERA. 

The formalism of the entanglement renormalization and the MERA were 
introduced in Refs. [5, 6]. Algorithms to approximate ground states have been 
described in Refs. [7, 8, 9]. Two dimensional systems have been explored in 
Refs. [10, 11, 12, 13, 14, 15, 16], including scalable simulations in interacting 
systems [13, 14] and analytical results for systems with topological order [15, 
16]. Scale invariant systems, including non-critical and critical fixed points of 
the RG flow, have been studied in Refs. [5, 6, 9, 10, 11, 15, 16, 17, 18, 19, 
20, 21]. Extensions to two-dimensional systems with fermionic and anyonic 
degrees of freedom have been proposed in Refs. [22, 24, 23, 25, 26]. 
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1.2 Coarse-graining and ground state entanglement 

Let us consider a system, such as a quantum spin model, that can be described 
by a lattice C in D spatial dimensions. For simplicity, most of the present 
derivations will involve a lattice in one spatial dimension. However, one of 
the highlights of entanglement renormalization is that it also applies to higher 
dimensional cases, to which we will refer occasionally throughout the chapter. 

The microscopic degrees of freedom are placed on the N sites of C, with 
each site being described by a vector space V of finite dimension d. The model 
is further characterized by a Hamiltonian H that decomposes as the sum of 
local terms, that is, of terms that act on a small number of neighboring sites. 
Except for the discussion on boundary critical phenomena in Sect. 1.5, we 
assume that the model is invariant under translations. 

Our ultimate goal is to be able to compute low energy properties of the 
system, which we will assume to be in the ground state |* G s) S W® N of H. 
Let 0\ , 02 , • • • ,Ok be arbitrary local operators acting on different parts of the 
lattice. Then wc would like to compute quantities such as 

(oi0 2 • • • o fc )^ GS = (* GS |oio 2 • • • o fc |* GS ), (1.1) 

since from these expected values one can predict how the system reacts to 
arbitrary external probes*. However, due to the exponential growth in N of 
the dimension of V® N t computing |^gs) by diagonalizing H is only affordable 
for very small systems. In order to address larger systems, we need a better 
plan. The strategy that we will pursue here is to build a transformation 
that removes short-distance degrees of freedom from the lattice model (C, H), 
which is mapped into an effective lattice model (£', H') such that 

(°i°2 • • • °fc)* GS = (oi02---o fe )* GS , (1.2) 

where \W GS ) is the ground state of H' and the operators o[, o' 2 , • • • ,o' k result 
from transforming o\ , 02 , • • • , Ok ■ 

Building an effective description can be helpful in several ways. On the one 
hand, in the case of a finite system, where the model (C,H') has a smaller 
Hilbert space dimension than the original model, diagonalizing H' is compu- 
tationally more affordable than diagonalizing H. As a result, larger models 
might be addressed. In conjunction with finite size scaling techniques, this 
already constitutes a powerful numerical route to study e.g. quantum critical 
phenomena. On the other hand, and closer to the goals of this chapter, a 
coarse-graining transformation that properly removes short-distance degrees 
of freedom could be used to investigate how H changes under scale trans- 
formations. Then, with the help of powerful RG ideas, we might be able to 
evaluate Eq. 1.1 directly in the thermodynamic limit. 



*We will only consider equal-time correlators. 
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FIGURE 1.1 

(i) Coarse-graining transformations characterized by an isometry w that maps 
blocks of three sites of lattice C into single sites of a coarse-grained lattice C . 

(ii) Graphical representation of w^w = hfi. (iii) An operator o supported on 
e.g. two blocks of sites of C becomes a two-site operator o', cf. Eq. 1.5. 

1.2.1 A real space coarse-graining transformation 

Following the seminal works of Migdal, Kadanoff and Wilson in real space 
RG [1, 2, 3], we will proceed by coarse-graining blocks of sites of C into single 
sites of £'. For concreteness, we divide C in blocks of three sites and, as 
in Wilson's numerical renormalization group (NRG) approach, we implement 
the coarse-graining by means of an isometry w, 

ui:V'^V® 3 , w^w = Iy>, ww< = P, P 2 = P, (1.3) 

where V® 3 is the vector space of the three sites, V is the vector space of the 
effective site, Iy is the identity operator in V and P is a projector onto the 
subspace of V® 3 that is preserved by the coarse-graining. Fig. 1.1 illustrates 
this transformation, which defines an effective lattice £ made of N' = N/3 
sites, as well as an effective Hamiltonian H' given by 

H' = W f HW, W = w® N/:i . (1.4) 

A local operator o with support on r blocks of C becomes a local operator o' 
supported on r sites of £' given by 

o -> d = (w^ r ) o , (1.5) 

where the isometries act on the relevant sites of C. Notice that the support 
of local operators may shrink under coarse-graining, but it never expands. 
In particular, if H can be expressed as a sum of interactions between pairs 
of nearest neighbor sites, then H' will also contain at most nearest neighbor 
interactions. 

The above transformation may not preserve the properties of the original 
ground state, in the sense of Eq. 1.2. As explained by White as part of 
his density matrix renormalization group (DMRG) [27], in order to preserve 
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ground state properties, the isometry w must retain the whole support of the 
ground state reduced density matrix p on the three-site block, 

P = tr B |* GS )(* GS |, (1.6) 

where E denotes all the sites of C not contained in the block. Let 

x _ 

P=2pa|*aX*a|, ^Pa = 1, Pi > P2 > ■ ■ ■ > 0, (1.7) 

a— 1 a 

be the eigenvalue decomposition of p. Then White's rule consists in choosing 
w such that the projector P in Eq. 1.3 corresponds to the support of p, 

x 

P = unJ = ^ l*aX*«l- White's rule (1.8) 

a=l 

In particular, the dimension of V' is the number \ of non-zero eigenvalues* of 
p. If Eq. 1.8 holds, then it can be seen that the ground state \W G3 ) of H' is 
given by |*' GS ) = W+|* os ), whereas WWt|* QS ) = |* os ). It follows that* 

(*gs|oio 2 • • • o fc |* GS ) = (y GS \(WW^)o l02 ■ ■ ■ o k (WW^ GS ) (1.9) 
= «* GS |W0 (^oioa • • • o k W) (W+|* os )) = (* GS K4 • • • O ' fe |* GS >,(1.10) 

which indeed amounts to Eq. 1.2. 

By iteration, from the original model (C^, H^) = (C,H) we can now 
build a sequence of increasingly coarse-grained lattice models 

(£(°) )jff (0)) ^ ^ (£( 2 ),tf( 2 )) (1.11) 

where lattice £( T+1 ) results from coarse-graining lattice C^ T \ such that ground 
state properties are preserved at each step, 

/ (r+l) (r+l) (r+l)> _ , (r) (r) (t)v , 

\°1 °2 °fe '*£s ~ \ 1 2 °fe 

How well does this coarse-graining scheme remove short-distance degrees of 
freedom? Is it viable in practice? To answer these questions, it is useful to con- 
sider the sequence of ground state reduced density matrices {g^ , g^ , g^ , • • • } 
corresponding to increasingly large blocks of sites of C, where g^ corresponds 
to 3 T sites. White's rule, which guarantees that ground state properties are 
preserved, also implies that the vector space dimension \^ of a site of lattice 
£( T ' must be given by the rank of g^ T \ Thus, our next step is to characterize 
the sequence of effective dimensions {x^°\ X^ 2 \ ' ' ' }■ 



tin practical calculations, one often neglects contributions from the eigenvectors \^ a ) with 
smallest weights p a , at the price of introducing small errors in Eq. 1.2. Most of the subse- 
quent discussions also apply in the case of approximate coarse-graining transformations, 
t For simplicity, we assume that the local operators o\ , o' 2 , • • ■ ,o' k arc supported on different 
sites of C , so that W^oi02 ■ ■ ■ o^W = o^o^ ■ ■ ■ oi. If e.g. operators Oj and o' 2 had overlap- 
ping support, then they would be fused into a single local operator o' = («>t® r ) o\02 (ui® r ) 
supported on some larger number r of sites. 
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1.2.2 Ground state entanglement 

The rank \ of the ground state reduced density matrix p of a block of sites 
depends on the amount of entanglement between that block of sites and the 
rest of the lattice, as can be seen from the so-called Schmidt decomposition, 
which expands the ground state |^gs) m terms of the eigenvectors \*5f a ) of p 
(Eq. 1.7) and some orthonormal set of states |$ a ) for the rest of the lattice, 

x 

|*GS> = ^V / ^l*a>® |*a>- (1-13) 
a=l 

When % = 1, corresponding to a reduced density matrix p = in a 

pure state, the ground state \^ G s) = l^i) ® |$i) factorizes into the product 
of individual states for the block and for the rest of the lattice. Instead, if 
X > 1, the block and the rest of the lattice are entangled. The amount of 
entanglement between these two parts can be measured with the entanglement 
entropy, namely the von Neumann entropy of the reduced density matrix p, 

x 

S(p) = -tr(plog(p)) = -J2Palog( Pa ). (1.14) 

a=l 

This measure vanishes for a product state and its maximum occurs for the flat 
probability distribution p a = 1/x, where S = —^2 a l/xlog(l/x) = l°g(x)- 
Therefore we always have x > exp(S'). In subsequent discussions we will 
assume for simplicity that 

X «exp(5). (1.15) 

In one spatial dimension, the entanglement entropy of a block of I contiguous 
sites typically increases with / until I becomes of the order of the correlation 
length £ in the system, at which point it saturates to some value S m&K , whereas 
it diverges logarithmically at a quantum critical point 

S(l) < 5 max , ID non-critical (1.16) 

5(0 « -log Z, ID critical (1.17) 

6 

where c is the central charge of the CFT that describe the infrared limit of 
the quantum phase transition. In two spatial dimensions, the entanglement 
entropy of a square block of I x I sites typically grows proportional to the size 
of the boundary of the block, 

S{l)Kicd, 2D (1.18) 

where a is some constant that depends on the model and grows with £, and 
where there might be logarithmic multiplicative corrections in some systems 
with quasi-long range correlations. More generally, in the ground state of a D 
dimensional model, the entanglement between a block of l D sites and the rest 
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of the system scales according to the so-called area law S(l) s=s al 0-1 , that 
is, as the size of the boundary of the block, with Eqs. 1.16 and 1.18 being 
particular cases of this expression. 

The area law translates into an approximate expression for the effective 
dimension x ■ Recall that one site of C^ T > accommodates I = 3 T sites of C 
(or I 2 = 3 2t sites if we had chosen to coarse-grain blocks of 3 x 3 sites in two 
spatial dimensions). Eqs. 1.16-1.18 combined with Eq. 1.15 lead to 

X (t) < Xmax ~ e 5max ID non-critical (1.19) 
x {t) w l c/6 w e T ID critical (1.20) 

x (t) » e al sa e eV¥ 2D (1.21) 

Now that we have an expression for the scaling of x^ T \ we can analyze both 
its origin and its implications. 



1.2.3 Accumulation of short-distance degrees of freedom 

From a conceptual viewpoint, the growth of x^ reveals an important flaw of 
the present coarse-graining scheme. Namely, that it fails to remove some of 
the short-distance degrees of freedom, which remain (and even accumulate) 
over successive iterations. 

To understand this point, let us assume that two contiguous sites r, s G C 
are in an entangled state, say (|l r l s )+|2 r 2 s ))/\/2. We compare two situations, 
depending on how the blocking in Fig. 1.1 affects these sites: 

(i) If both sites belong to the same block, then they do not contribute to 
the entanglement between the block and the rest of the lattice, or what is the 
same, to the spectrum of the reduced density matrix p of the block, and are 
therefore removed during the coarse-graining. 

(ii) If the same two sites belong to two adjacent blocks, then they are 
entangled across the boundary of the blocks, contribute to the spectrum of p 
for each block, and thus will be preserved by the coarse-graining, see Fig. 1.3. 

In other words, short-distance degrees of freedom, entangled at scales smaller 
than the size of a block, arc treated differently depending on whether they are 
entangled within a block or across the boundary between blocks. In the last 
case, the coarse-graining transformation fails to remove them. As a result, the 
approach cannot generate a proper RG flow. For instance, two Hamiltonians 
Hi and H2 describing the same phase but differing in short-distance details 
will remain different under successive iterations of the coarse-graining: some 
of the short-distance details that distinguish the two models will persist, and 
two different fixed points of the RG flow will be reached. 

In spite of this shortcoming, is the scheme useful for practical computations? 
The answer strongly depends on the class of problems under consideration. 
It is a definitive 'yes' for systems in one spatial dimension. White's DMRG 
[27] (to which this scheme is closely related) revolutionized the computational 
scene in the early 90's by providing numerical access, with unprecedented 
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FIGURE 1.2 

(i) Coarse-graining transformation with disentanglers u, acting on two con- 
tiguous sites of C across the boundary of two blocks, followed by isometrics 
w, which map a block of three sites into a single site of the coarse-grained 
lattice C . (ii) Graphical representation of uv) = 7y®2. (iii) A local operator 
o supported on two contiguous sites of C is mapped into an effective operator 
d on two contiguous sites of £' . The linear transformation o — > d = A(o) 
is described by a superoperator (i.e., a linear map in the space of operators) 
referred to as the ascending superoperator A. 



accuracy, to ground state energies, local observables and correlators. The 
immense success of DMRG can be understood from the saturation of 
for non-critical systems, Eq. 1.19, since the computational cost is a small 
power of x • Even critical systems can be studied, in spite of Eq. 1.20, by 
considering reasonably large chains, of up to thousands of sites, and applying 
finite-size scaling techniques. Instead, the situation is much grimmer in two 
spatial dimensions. There the double exponential growth of x^ T ' in Eq. 1.21 
implies that only small systems can be addressed (see e.g. the discussion in 
Ref. [28]) and there is a real need for alternative approaches. 



1.3 Entanglement Renormalization 

Entanglement renormalization was proposed as a means to consistently re- 
move all short-distance degrees of freedom during coarse-graining. It leads to 
a significant reduction in computational costs, making scalable calculations 
feasible even in two dimensions. More important to us, it produces a RG flow 
with the expected properties. 
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1.3.1 Disentanglers 

Let us modify the previous coarse-graining scheme as follows. First we apply 
a unitary transformation u, a disentangle^ on pairs of contiguous sites of C, 

u:V® 2 ^V® 2 , u ] u = uu ] =7 V ®2, (1.22) 

where each site belongs to one of two adjacent blocks, so that u acts across the 
boundary of the two blocks, see Fig. 1.2. Then, as before, we use isometrics 
w to map the blocks of 3 sites of C into single sites of the lattice Let 
U = u® N / 3 be a unitary transformation that contains disentanglers u acting 
on the boundaries of all the blocks. Then the effective Hamiltonian in £ reads 

H' =(W^U ] )H{UW). (1.23) 

In order to guarantee the preservation of ground state properties, the isome- 
tries w are chosen according to White's rule, namely such that they retain the 
support of the reduced density matrix p on the block, 

P = tr £/6iocfe (C/t|* GS )(* GS |f/). (1.24) 

Notice, however, that \^ GS ) has been transformed according to the disen- 
tanglers before the density matrix p is computed. By properly choosing the 
disentanglers u, the modified ground state |^ GS ) = l^|\1/ gs ) will have less en- 
tanglement than l^cs)- Specifically, the disentanglers are capable of removing 
short-range entanglement across the boundary of the blocks, see Fig. 1.3. 

By iterating the transformation we can once more build a sequence of in- 
creasingly coarse-grained models (Fig. 1.4) 

(£<°),ij(°)) — ► (£«, H^) — ► 2 \H^) (1.25) 

As evidenced by abundant numerical and analytical results (Refs. [5] -[26]), an 
important consequence of the use of disentanglers is that the dimension of 
the effective sites no longer needs to increase with r in order to (approximately, 
but accurately) fulfill White's rule. This is the case in one spatial dimension, 
both for critical and non-critical systems, cf. Eqs. 1.16-1.17, as well as in 
two spatial dimensions, provided the area law of Eq. 1.18 does not have 
logarithmic corrections. 

1.3.2 Ascending and descending superoperators 

An operator o supported on two contiguous sites of C is now mapped into an 
operator d supported on two contiguous sites of £ by means of the ascending 
superoperator A (Fig. 1.2), 

o^o' = ((t«t® 2 )( u t)) o (u{w® 2 )) = A(p). (1.26) 



§The distinction between disentanglers u and isometrics w is useful for pedagogical purposes. 
However, more generally one can consider isometric transformations v : V®" 1 i— ► V®" 2 , with 
v^v = -Tyigni , ni < ri2, that both disentangle the ground state and coarse-grain the lattice. 
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FIGURE 1.3 

To illustrate the role of disentanglers u, a simple example is considered where 
the ground state \^> GS ) of a one-dimensional lattice C factorizes into the prod- 
uct of entangled states (|l r l s ) + |2 r 2 s ))/v / 2 involving only two nearest neigh- 
bor sites r,s € C (on both sides of the boundaries between blocks) and states 
of single sites (in the interior of each block), (i) When only isometries are 
used, entanglement across the boundary of the blocks is preserved in £' . (ii) 
By using disentanglers u such that transform the state (|l r l s ) -I- |2 r 2 s ))/v2 
of the two boundary sites r,s £ £ into an unentangled state, e.g. |l r l s ), 
entanglement across the boundary of the blocks can be removed before the 
isometries are applied, and the ground state of £ has no entanglement. 




FIGURE 1.4 

Diagramatic representation of two iterations of the coarse-graining transfor- 
mation, producing a sequence of increasingly coarse-grained lattices , £W , 
and CP^ 1 . At each iteration is chosen as to remove short-range entangle- 
ment and u/ T ) follows White's rule. The ascending superoperator A maps 
o( T ) into o( r+1 ) while the descending superoperator V maps p^ T+1 ^ into p^ T \ 
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Notice that locality is again preserved under the coarse-graining transfor- 
mation, but in this case it results from a compromise between two opposing 
forces: disentanglers u expand the support of local operators, while isometries 
w compress it. The balance corresponds to two-site supports, to which both 
smaller and larger supports tend under coarse-graining. This is why two-site 
supports play a dominant role in the present discussions. In the case of a 
sequence of coarse-graining transformations (Fig. 1.4) we can investigate how 
the two-site operator = o changes as we increase the scale of observation: 

o (0)^ o (D^ o (2)_>... (1.27) 

On the other hand, the descending superoperator V, 

P^P = tr SlS2 (K 2 )«) d , (1.28) 

where p and p' are supported on the same two sites as o and d in Eq. 1.26 
and the partial trace eliminates two superfluous sites s\,si G C, produces a 
fine-grained reduced density matrix p from a coarse-grained reduced density 
matrix p' , such that 

tr(po) = tr(p'o'), (1.29) 

for all possible two-site operators oT By iteration, the descending superop- 
erator allows us to obtain a sequence of two-site density matrices 

> ,(2)2^(1)2^(0) (L30) 

flowing from coarser to finer lattices and therefore monitor how the ground 
state two-site density matrix changes with the observation scale. 

1.3.3 Multi-scale Entanglement Renormalization Ansatz 

The multi-scale entanglement renormalization ansatz (MERA) is a variational 
ansatz for pure states |^) € V® N of C that results in a natural way from the 
above coarse-graining transformation. It consists of a tensor network that 
collects all the disentanglers and isometries used to produce the sequence of 
coarse-grained lattices in Eq. 1.25, see Fig. 1.5. 

In the case of a translation invariant system, we need to specify one dis- 
entangler and one isometry per layer, and these depend on 0(x 9 ) 
parameters, where q — 4 in the one-dimensional scheme discussed here. Since 
there are O (log (TV)) layers of tensors, the MERA depends on 0(x q log(iV)) 
parameters, a result valid also in two spatial dimensions for an appropriate 
value of q. 



If In other words, the descending superoperator V is dual to the ascending superoperator 
A, tv(p'A(o)) = tr(V(p')o) for all p' and o. 
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FIGURE 1.5 

Multi-scale entanglement renormalization ansatz made of disentanglers and 
isometries corresponding to two iterations of the coarse-graining transforma- 
tion in Fig. 1.2. Notice the periodic boundary conditions. 



Therefore, the MERA is an efficient representation of certain states of lat- 
tice C. The interest in this ansatz resides in that abundant numerical and 
analytical results demonstrate that it can be used to accurately approximate 
the ground state I'J'gs) of a local Hamiltonian H in large (and even infinite) 
lattice systems in one and two spatial dimensions. Even at criticality the er- 
ror in ground state expected values, Eq. 1.1, decays exponentially with the 
refinement parameter \- While it is still unclear exactly what Hamiltonians 
have ground states that can be well approximated, a MERA with a fixed 
value of x naturally reproduces a scaling of entanglement compatible with 
Eqs. 1.16-1.18, suggesting that it might be able to accurately approximate 
ground states whose block entanglement fulfills the area law. 



Given the ground state I'J'gs) of a local Hamiltonian H, we have pro- 
vided an intuitive description of the role played by the tensors that form the 
MERA. Namely, a disentangler u 'removes short-range entanglement from 
I'J'gs)', while an isometry w 'coarse-grains a block of sites so as to preserve 
the support of the ground state reduced density matrix'. These descriptions 
could perhaps allow us to compute all it's and w's of the MERA for I'J'gs) if we 
already knew I'J'gs)- Algorithmically, of course, it is absurd to assume that we 
already know I'J'gs) if our goal is to compute I'J'gs)- In practical MERA cal- 
culations, an approximation to I'J'gs) is obtained by minimizing the expected 
value ('J'|i7|'J'), where |\&) is constrained to be a MERA. This is done starting 
with random u's and w's and iteratively optimizing them so as to reduce the 
energy, with a cost that scales as 0(x q l°g(A0) P er optimization step (q' = 7 
in the present scheme). We refer to the literature for further details. 
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1.4 The Renormalization Group picture 

From now on we consider an infinite lattice C and a Hamiltonian H that is 
the sum of nearest neighbor terms 

H = ^2h(s,s + l), h{s,s + l) = h, (1.31) 

s 

where the two-site operator h is the same on all pairs of nearest neighbor 
sites, and thus h completely characterizes H. We assume that disentanglcrs 
and isometries have been properly optimized, e.g. by minimizing the expected 
value (^\H\^), and concentrate on analyzing the properties of the resulting 
coarse-graining from the perspective of the renormalization group. 

1.4.1 A real space RG map 

Let us introduce the average ascending superoperator A h , 

> = A H0)= AhL{0)+Ahc{0)+AU °\ (1-32) 

3 

that averages over the three inequivalent ways (left, center and right) in which 
o can be coarse-grained, and where the superscript h highlights the depen- 
dence on the Hamiltonian H — recall that the ascending superoperator is built 
with the disentanglers and isometries corresponding to the ground state \$> GS ) 
of H . Under coarse-graining, H is mapped into an effective Hamiltonian H', 

H' = 3^2h\s,8 + l), h'(s,s + l) = h', (1.33) 

s 

where the constant two-site operator h! is obtained from h by the RG map 1Z, 

ti = K(h)=A h (h). (1.34) 

Notice that, by construction, IZ(h) is a non-linear function of h. 

Given the two-site interaction term = h of the initial Hamiltonian 
we can now build a sequence of two-site interactions 

h<°> A h {1) A h {2) A • • • (1.35) 

where characterizes the Hamiltonian = 3 r ^ s (s, s + 1) of the 
coarse-grained lattice £( T \ This defines a discrete RG flow in the (x 4 - 
dimcnsional) space of possible two-site interactions that we can use to study 
how the Hamiltonian of the system changes with the observation scale. 
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FIGURE 1.6 

The RG map 7Z(h) defines a RG flow in the space of two-site interactions h. 
The scaling superoperator <S(o), Eq. 1.41, transforms two-site operators o at 
the fixed points of this RG flow. 



1.4.2 Properties of the RG map 

Let us briefly discuss a few properties of the RG map IZ(h): 

(i) Proper RG flow. — Abundant numerical evidence suggests that properly 
chosen disentanglers indeed succeed at removing all short-distance degrees 
of freedom. As a result, if two Hamiltonians Hi and H 2 correspond to the 
same phase, multiple applications of the RG map take and into the 
same fixed-point interaction h* , that is h\ — h\ (up to trivial changes of local 
basis). 

(ii) Nearest neighbor interactions. — As mentioned earlier, the present coarse- 
graining transformation does not generate long-range interactions starting 
from a short-ranged Hamiltonian (in contrast e.g. with momentum-space RG 
methods). In addition, a Hamiltonian containing arbitrary short-range inter- 
actions (well beyond nearest neighbors) can be reduced, after a few iterations, 
to a Hamiltonian with only nearest neighbor interactions". Therefore the two- 
site RG map IZ(h) can be used to study arbitrary phases with short-range 
interactions. In two spatial dimensions, the analogous construction leads to 
an RG flow for four-site Hamiltonians. 

(Hi) Unbiased RG map. — Notice that the space V for an effective site of £ 
is not chosen a priori on the basis of heuristic arguments, nor is some specific 



II A coarse-graining step with trivial disentanglers is required (once) in order to eliminate 
some of the next-to-nearest neighbor interactions. 
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form of the effective Hamiltonian H' (with a few free parameters to be fixed) 
imposed. The approach simply searches, through an energy minimization, the 
subspace V® N / 3 of the total Hilbert space Y® N that best approximates the 
low energy subspace of H. The only hypothesis that is made concerns the 
structure of the subspace Y® N ' 3 , assumed to factorize as a tensor product 
of spaces V that are related to the original factorization Y® N through the 
discntanglers u and isometries w. Numerical evidence (provided e.g. by the 
study of critical phenomena in the next section) suggest that the low energy 
subspace of many local Hamiltonians of interest indeed has this structure. 

1.4.3 Fixed points of entanglement renormalization 

Following the RG program, our next step is to study the fixed points of the 
flow generated by the RG map 1Z, corresponding to Hamiltonians that do not 
change when we modify the scale of observation. That is, we are interested 
in models with two-site interactions h* such that 

1Z(h*) = Xh* scale invariance (1.36) 

where A = 3~ Ah * is a proportionality constant to be discussed later on. 

A fixed point is characterized by one discntangler u and one isometry w 
that are repeated at all length scales, i.e. — u and = w for all r. 
Thus, the scale invariant MERA, specified by this pair (u,w), depends on 
just 0(x 4 ) parameters and offers an extremely compact description of fixed- 
point ground states of an infinite lattice. It can be obtained with an algorithm 
whose cost formally scales as 0(x 7 ) [18]. 

Since a finite correlation length £ shrinks under coarse-graining, scale in- 
variance requires that £ = or £ = oo, corresponding to non-critical and 
critical fixed points. The present approach leads to a natural characterization 
of RG fixed points in terms of the entanglement of their ground states and 
the MERA representation: 

(i) Non-critical fixed-point ground states (£ = 0) may be unentangled or en- 
tangled. Symmetry-breaking phases usually have an unentangled (or product) 
fixed-point ground state, represented with a trivial scale invariant MERA with 
X = 1, whereas the fixed-point ground state in topologically ordered phases 
is entangled. The results of [15, 16] strongly suggest that a scale invariant 
MERA with finite x can exactly represent such entangled ground states. 

(ii) Critical fixed-point ground states (£ = oo) are always highly entangled. 
An exact MERA representation requires an infinite dimension \, as will be 
discussed in Sect. 1.5, but we will also see that accurate estimates of critical 
properties can be obtained from a scale invariant MERA with finite x- 

In order to further characterize a fixed point, one usually linearizes 1Z at 
h = h* by considering its derivative with respect to a small perturbation eo, 

,. H(h* + eo) - TZ(h*) ,. A h * +eo (h* + eo) - A h * (h*) 

hm = hm . (1.37) 
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We conjecture** that this derivative is given by A h (o), which dictates how 
two-site operators o transform under coarse-graining at the fixed point h = h* , 
and which we will call scaling superoperator S, 

S(o)=A h *(o). (1.38) 

Recall that one can distinguish between stable fixed points, where any per- 
turbed interaction h* + to will flow back to h*\ and unstable fixed points, 
where relevant perturbations o exists such that h* + to will flow away from 
h*. It may also be worth recalling that the notion of stability is relative to 
what perturbations o are available, and thus may depend on symmetry con- 
siderations. As illustrated in the next section for critical fixed points, from 
the scaling superoperator 5 one can identify all relevant perturbations of a 
fixed point and therefore analyze its stability. 



1.5 Quantum criticality 

In this section we apply the scale invariant MERA to the study of critical 
fixed points of the RG flow, that is to quantum critical points. In this way 
we describe an important application of entanglement renormalization, while 
also demonstrating the validity of the approach. For concreteness, we analyze 
the quantum Ising chain with critical transverse magnetic field, 

s s 

whose well-known critical properties are described by a (1+1) conformal field 
theory (CFT) [29]. 

1.5.1 Scaling operators and critical exponents 

Most properties of a quantum critical point can be extracted from the spectral 
decomposition of the (linearized) RG map at the corresponding fixed point of 
the RG flow, see e.g. [4]. In the present case, the spectral decomposition of 
the scaling superoperator S [17, 18], 

S(») = ^2 <Watr(<£ a »), tT(4> a <f>/}) = 5 a p, (1.40) 

a 

readily provides us with the scaling operators of the lattice model, namely 
those operators <f> a on two contiguous sites that transform into themselves 

"This conjecture would be true if \im t ^(A h + eo (/i*) — A h (h*))/e=0, which at a critical 
point is both plausible and compatible with numerical tests. 
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FIGURE 1.7 

Scaling dimensions of the quantum critical Ising model extracted from the 
bulk and boundary scaling superoperators of a \ — 16 MERA. 



FIGURE 1.7 

Scaling dimensions of the critical quantum Ising chain. 
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under the coarse-graining transformation, 

S{dp a ) = KK, A Q =3- A ", (1.41) 

as well as their scaling dimensions A Q = — log 3 (A Q ). An example of scaling 
operator is the identity operator fa = I, 

S(l) = I, Ai = 0, (1.42) 

as can be seen by using u*u = I and w^w = I in Eq. 1.26. Its dual operator 
4>i in Eq. 1.40 is the ground state reduced density matrix fa = p for a block 
of two contiguous sites. A second example of scaling operator is of course 
the two-site fixed-point interaction h* , see Eq. 1.36, with scaling dimension 
A h * =2. 

Fig. 1.7 shows the 12 smallest scaling dimensions for the Ising model, 
obtained with a scale invariant MERA with \ = 16. As expected, they 
appear organized in the conformal towers of three scaling operators known 
as the identity I, the spin a, and the energy e, which correspond to the 
primary fields of the Ising CFT [29]. Recall that all critical exponents of 
a model can be obtained from the scaling dimensions of its primary fields. 
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FIGURE 1.8 

One-point, two-point and three-point correlators for one-site scaling operators. 
After \og 3 (r xy ) iterations of the coarse-graining transformation, two operators 
separated by r xy sites become nearest neighbors. Hence the polynomial decay 
of correlations. 



For instance, for the Ising model the exponents v and r\ are v = 2lS. a and 
f] = 2 J A , and the scaling laws express the critical exponents a,(3,j,S in 
terms of v and r\. We can also identify e.g. all the relevant perturbations of 
the model, corresponding to scaling dimension A < 2, as well as its marginal 
perturbations, with A = 2, which are listed in Table 1.1. 



1.5.2 Correlators and the operator product expansion 

Extracting correlators from a scale invariant MERA is also relatively easy. 
In terms of the one-site scaling operators ipa (see [18] for details), one-point, 
two-point and three-point correlators read, 



(ip a (x)ip (y)) 



Cq/3 



Pi 



A Q +A,3-A 7 



Ah+A^-Ao 



A^+A 



(1.43) 



~^a7' (I- 44 ) 



where the sites x, y, z have been chosen conveniently and r xy = \x — y\, r yz = 
\y - z\ and |. The coefficients C a , C a p and C a ^ are given by 

C a =tr(4> a p) = 5 a i, C a p = tr ((ijj a ® ipp)p) , (1.45) 

C aM = 2 a "+ a ^- a Ht ((V» q ® V/3 ® ^)p) , (1.46) 
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where p simultaneously denotes the reduced density matrix on one, two and 
three sites respectively. Thus, as expected in a critical system, the scale invari- 
ant MERA produces polynomial correlators. How this occurs is very intuitive. 
Consider for instance the two-point correlator (ip a (x)il>f}(y)) for r xy — 3 r . As 
illustrated in Fig. 1.8, operators ip a and ipp become nearest neighbors af- 
ter r = \og 3 (r xy ) iterations of the coarse-graining transformation. Since this 
transformation maps ip a into tp' a = 3~ Aa ip a , each iteration contributes a fac- 
tor 3- A £»- A /3 to the correlator, with 

^-Ac-Apy _ ^-Aa-A^loga^y) _ 3^3(7-^ " Aq ~ A £ ) _ ]_ ^ (1.47) 

Txy a 13 

which explains its scaling. Finally, in order for the correlator (ip a (x)ipf}(y)) 
to be non-zero, the tensor product of the scaling operators ip a (g) ipp must fuse 
into the two-site identity operator I, which occurs with amplitude C a p. 

The coefficients of Eqs. 1.45-1.46 are analogous to those that appear in 
CFT. A proper choice of normalization of the fields leads to C a p — 5 a p^ 
whereas coefficients C a p 1 define the operator product expansion (OPE), which 
for the primary fields of the Ising CFT reads 

/1CFT c /^CFT ^ pCFT /^CFT /^CFT ft /-I AQ\ 

A MERA with x = 16 reproduces all these OPE coefficients with errors 
smaller than 3 x 10~ 4 . 

The above comparisons with exact results from CFT allows us to conclude 
that entanglement renormalization produces accurate estimates of the critical 
properties of a system. Once this has been established, the present approach 
can be used to actually identify which CFT describes a given quantum critical 
point by estimating the conformal data (central charge, conformal dimensions 
and OPE of primary fields) that completely characterize it [18]. The present 
discussion also indicates that the scale invariant MERA can be regarded as 
approximately realizing an infinite dimensional representation of the Virasoro 
algebra [29]. The finite value of \ implies that only a finite number of the 
quasi-primary fields of the theory can be included in the description. Fields 
with small scaling dimension, such as primary fields, are retained foremost. 
However, an exact description requires x to be infinite. 



1.5.3 Boundary critical phenomena 

So far we have restricted our attention to translation invariant systems. Crit- 
ical systems with a boundary can also be described by adding a boundary 
to the scale invariant MERA [21]. This is done by introducing a boundary 
isometry ui surf - at the boundary, while the pair (u, w) corresponding to the 
translation invariant case still represents the bulk, see Fig. 1.9. 
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FIGURE 1.9 

Scale invariant MERA with a boundary made of copies of the same boundary 
isometry w Buri and the tensors (u, w) of the bulk. 

FIGURE 1.9 

Scaling dimensions for free and fixed boundary conditions 
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The presence of a boundary modifies the bulk correlators of Eqs. 1.43-1.44. 
For instance, the expected value (tp a (x)} becomes non-trivial 

(^WU^, (1-49) 

where x is the distance to the boundary (in number of sites) and Cg™ f is some 
constant. Again, this result is very intuitive. Consider x = (3 T+1 — l)/2. 
As illustrated in Fig. 1.9, after r = log 3 ( 2x g" x ) rj log 3 (x) iterations of the 
coarse-graining transformation the bulk scaling operator ^/>^ ulk has reached 
the boundary, into which it fuses with amplitude Cq™ 1 ' ■ In addition, by di- 
agonalizing the boundary scaling superoperator S eutt that maps the boundary 
into itself at different length scales, we can extract the boundary scaling op- 
erators and their scaling dimensions, as well as the fusion rules between 
bulk and boundary operators. Fig. 1.7 and Table 1.2 show the boundary 
scaling dimensions for free and fixed boundary conditions. One can see that 
the number of conformal towers is smaller than in the bulk: for free boundary 
conditions, the primary fields left are the identity I and the spin a, whereas 
for fixed boundary conditions only the identity I remains. Once more, the 
results compare well with the exact solution provided by (boundary) CFT 
[29]. 
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1.6 Summary & Outlook 

In this chapter we have reviewed, step by step, the construction of a coarse- 
graining transformation for quantum many-body systems on a lattice that 
fulfills two natural requirements: it disposes of all short-distance degrees of 
freedom while it retains the properties of the ground state of a local Hamil- 
tonian. This is achieved by using disentanglers u, that remove short-range 
entanglement across the boundary of a block of sites, and isometries w, that 
coarse-grain these blocks of sites according to White's rule and thus preserve 
the support of the ground state reduced density matrix. 

This transformation can be used to compute ground state expected values 
of the form (oio 2 ■ ■ ■ o^) both in finite and infinite systems. In a finite system, 
we simply coarse-grain the lattice until the effective Hamiltonian can be nu- 
merically diagonalized. In an infinite system, we iterate the coarse-graining 
transformation until we become sufficiently close to a fixed point of the RG 
flow — that is, until we have eliminated all irrelevant perturbations in the 
original Hamiltonian. Then we expand local operators in terms of the scaling 
operators of the theory, whose expected values can also be determined. In 
a symmetry-breaking phase, the fixed point corresponds to an unentangled 
ground state, and local operators can be trivially evaluated. In a critical 
phase, the fixed point corresponds to a highly entangled ground state, but a 
scale invariant MERA can still be used to approximately compute the scaling 
operators and evaluate their expected values. 

Notice that the scale invariant MERA allows us to make precise the notion 
of scale invariance in a lattice system. In order to get rid of the dependence on 
the original lattice spacing, one typically invokes the continuum limit. Here, 
instead, we have considered a coarse-graining such that all short-distance 
degrees of freedom are consistently discarded. The ground state of the lattice 
is then invariant under changes of scale if its effective description is locally 
identical to the original one. We envisage that this lattice version of scale 
invariance will become a useful testing ground for ideas and problems that 
are harder to analyze in the continuum limit. Recent work on the holographic 
principle might be a first step in this direction [30] . 

We conclude this chapter with a suggestive interpretation of the MERA. 
Notice that an isometry w : V — > V® 3 can be regarded as a unitary transfor- 
mation v : V® 3 — > V® on three sites where two of these sites are initially in 
some fixed state |0), that is to = v(\0) <g> |0)). The MERA can thus also be 
understood as a quantum circuit, made of gates v and u, that implements a 
unitary transformation V such that an unentangled state (0)®^ becomes the 
ground state \\& G s) of a local Hamiltonian H, see Fig. 1.10. Thus, ground 
states of local Hamiltonians can be regarded as the result of a quantum com- 
putation that flows from the largest available length scales (i.e., the size of 
the system) to the smallest length scales (distance between sites of C). At 
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FIGURE 1.10 

The MERA seen as a unitary transformation V such that |^ G s) = V\0)® N . 

each length scale, unentangled wires in state |0) are added and subsequently 
entangled with the rest of the system by means of the gates v and u. 

The quantum circuit interpretation is quite insightful. For instance, it pro- 
vides an explicit recipe to experimentally prepare the state I^gs) by entangling 
quantum systems, and it has been instrumental in generalizing the MERA to 
deal with fermionic and anyonic degrees of freedom [22, 24, 23, 25, 26]. It 
also highlights the reversible character of the coarse-graining transformation. 
During coarse-graining, short-distance degrees of freedom are not actually 
'thrown away', but rather 'stored' in disentanglers and isometries, and can 
be re-incorporated in the picture later on. Indeed, as we have seen, from the 
reduced density matrix describing the state of the system at a given length 
scale, we can recover the state at a smaller length scale by using the de- 
scending superoperator T>. Finally, this viewpoint emphasizes the astonishing 
structural similarities shared by most ground states of local Hamiltonians. 
Consider two systems in two completely different symmetry-breaking phases 
but with the same correlation length £. Both ground states become essentially 
disentangled after 0(log(£)) iterations of the coarse-graining. Therefore, the 
two ground states can be prepared from the same unentangled state by using 
the same quantum circuit, made of 0(log(£)) rows of identically wired gates, 
just by modifying the specific gates v and u applied in each case. 
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